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■ Abstract 

. We consider tiie feedback stabilization of a variable profile for an ensemble of non interacting half spins described by the 

■ Bloch equations. We propose an explicit feedback law that stabilizes asymptotically the system around a given arbitrary target 

■ profile. The convergence proof is done when the target profile is entirely in the south hemisphere or in the north hemisphere of 
I the Bloch sphere. The convergence holds for initial conditions in a neighborhood of this target profile. This convergence is 
, shown for the weak topology. The proof relies on an adaptation of the LaSalle invariance principle to infinite dimensional 

■ systems. Numerical simulations illustrate the efficiency of these feedback laws, even for initial conditions far from the target 

■ profile. 
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■ 1 Introduction 

Ensem ble controllability as introduced in lLi and Khaneial 
! (I2009D is an interesting control theoretic notion well 
' adapted to nuclear magnetic resonance (NMR) sys- 
] tems (see, e.e;..[Li and Khan eja (2006) and the reference 

■ herein) . In iBeauchard et al.l (|2010() some controllability 

■ issues of such NMR systems are investi gated using open- 
loop c ontrols involving Dirac-combs. In lBeauchard et al] 
(|201l[ ) such open-loop Dirac-combs are combined with 
Lyapunov stabilizing feedback to ensure closed-loop 
convergence towards a target profile that is one of the 
two steady-states, the south and north poles of the 
Bloch sphere. In this note, we extend this Lyapunov 
design to arbitrary target profiles and prove its local 
convergence for weak topology when the target pro- 
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file lies entirely in the south hemisphere or in the north 
hemisphere. 

We consider an ensemble of non interacting half-spins 
in a static field (0, 0, BqY in M.^, subject to a transverse 
radio frequency field {ui{t),U2{t),0Y in (the control 
input). The ensemble of half-spins is described by the 
magnetization vector M G depending on time t but 
also on the Larmor frequency uj = — 7-Bo (t is the gyro- 
magnetic ratio). It obeys to the Bloch equation: 

dM 

— (<,w) = {ui{t)ei + U2{t)e2 + ujes) A M{t,uj), (1) 
ot 

where — oo < < w* < -l-oo, w G (0;*,^*), (61,62,63) 
is the canonical basis of M^, A denotes the wedge product 
on R^. The equation (1) is an infinite dimensional bilin- 
ear control system. The state is the w-profile M, where, 
for every cj G {uj^,uj*), M{t,uj) G (the unit sphere of 
M.^). The two control inputs ui and U2 are real valued. 

We propose here a first answer to the local stabilization 
of an arbitrary profile: given an arbitrary target pro- 
file Mf : (u!^,,u!*) — 7- define an explicit control law 
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{ui{t, M), U2{t, M)), a neighborhood U of Mf (in some 
space of functions to be determined), a diverging se- 
quence of times (i„)„gN, such that, for every initial con- 
dition Af" e U, the solution of the closed loop system is 
uniquely defined and satisfies 



lim \\M{tr. 

i—>--\-oo 



Mf{. 



I {uj^ ■,^*) 



0. 



In this note, the Lyapunov feedback proposed in lBeauchard 
(|201l[ ) is adapted to provide a constructive answer to 
this question. Section 2 is devoted to control design and 
closed-loop simulations. In section 3 we state and prove 
the main convergence result, theorem 1. 



2 Lyapunov approach 

2. 1 Some preliminaries 

Let us recall the concept of a solution for (1) when 
the control input u contains Dirac distributions. When 
Ui,M2 € L]^^{M), then, for every initial condition Mq G 
w*), K^), the equation (1) has a unique weak so- 
lution M G C°{[Q,+oo),L^{{uj-^,w*),B?)). Denote by 
5{t — a) the Dirac distribution located at t — a. When 

ui = a6{t — a) + u\ and U2 — u\ where G L]g^{M), 
a > and a G (0, +cx3), then the solution is the classical 
solution on [0, a) and (a, -l-oo), it is discontinuous at the 
time t — a, with an explicit discontinuity given by an 
instantaneous rotation of angle a around the axis Rei 



of variable 

Mi{t,uj) := P{t)M{t,uj) where P{t) 
one gets the following dynamics 



f 1 \ 

e{t) 
\0 e{t)) 



(3) 



1 



\ 



cos(a) — sin(a) 
y sin(a) cos(a) J 



M{a 



The symbol ||.|| (resp. (., .)) denotes the Euclidian norm 
(resp. scalar product) on and the associated operator 
norm on 



dt 



{t,uj) = [ul{t)el+U2{t)e2 + e{t)uJ\^Ml{t,uJ). (4) 



The application of impulses at i = fcT, by changing the 
sense of rotation of the null input solution, is expected 
to reduce the dispersion in the closed loop system. Since 
M(i, w) = Mi(t,w) for every t G [2kT, {2k + 1)T], any 
convergence result on Mi{t) when t — > +00 provides a 
convergence result on M. 

The first step of the control design consists in putting 
the system (4) in driftless form. The new function 

M2(t,w) := eTip[a{t)ujS]Mi{t,uj) 



where 



a{t) :— / e{s)ds, S := 
Jo 



f io\ 

-10 




(5) 



solves 

dM2, 
~dt 



{t,uj) =^M,(t)[exp(a(t)wS')eJ AM2(t,w). (6) 



Since a{2kT) = 0,Vfc G N, any convergence on M2{t) 
when t — > +CXD provides a convergence on Mi {2kT) when 
k — !> +00. 

2.3 Transformation of the target profile 



2.2 Transformation into a driftless system 

As in lBeauchard et all ()2011[ ) we consider a control with 
an "impulse-train" structure 

+00 

ui ^ ui + "^Tr S{t-kT), U2 = {-iy^*'^U2 (2) 
fc=i 

where e{t) :— E{t/T), for some period T > and i?(7) 
denotes the integer part of the real number 7. The new 
controls ui , U2 belong to L]^^iM) . Considering the change 



The second step of the control design consists in trans- 
forming a convergence to a variable profile Mf into a 
convergence to the constant profile — 63, for which we 
developed tools in the previous work iBeauchard et all 
(|201l[ ). It relies on the following proposition. 

Proposition 1 There exists C > such that, for all 

Mf G iJi(K,'^*),S^); there existsR G H^{{uj^,u*), SO^iM.)) 

satisfying 

R{uj)Mf{uj) = -es, VwGK,a;*], (7) 
\\R\\h. ^C\\Mf\\H^. (8) 
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Proof: Let Mj e iJi((a;*, w*), S^) and set /(w) 
Mj:{uj) A Mf{uj). Denote by A{uj) the skew-symmetric 
operator defined by ^ M >^ f(uj) A M e Con- 
sider the Cauchy problem 



R = RA{uj) on [uj^,uj*] with R{uj^) = R^ 



where i?* is any rotation sending M/(cli,) to —63: 
RifMf{u)^:) = —63. Since H> is the solution i? 
is well defined, unique and belongs to i7^((a;,,w*),S^). 
Direct computations show that -^{RMf) — 0. Thus 
R{uj)Mf{uj) = -63. Moreover, II = 1 and 
||i?'(a;)|| = \\A{uj)\\ = ||M}(a;)|| for all uj e [uj.,uj*], 
which proves (8). □. 

Let us consider a target profile M/ e iJ^((w*,a;*),S^). 
Take R G cj*), S'03(R)) given by the above 

proposition. To any solution M2 of (6), we associate the 
function 

N{t,Lj) ;= i?(w)M2(t,w),Vw e [-a;,,cj*]. (9) 
This function solves the equation 

dN , 



dt 



i=l 



F(i,cj)e, AN{t,uj) (10) 



where 



F{t,Lu) := i?(w)exp(cr(t)w5). 



(11) 



The convergence of N(t, lu) to —63 ast +00 is equiva- 
lent to the convergence of A/2 toMf{uj) asi — ;> -foo. 



2.4 Lyapunov feedback 

Let us consider the following Lyapunov-like functional 



£(7V) := 



11^^ + ^311?^ 


1 


2 










dN 




duo 



+ l + {N,e-i))duj. (12) 



The function C is defined for any TV e ^((w*, w*), S^) 
and takes its minimal value on this space at the point 
N — —63 with £(—63) = 0. For any solution of (10), 
some computations show that 



dN , , fdF ^ , \ A 
-(.),(-(.,.)e.)ATV(.)) 

+ ^e3,(i^(t,w)ei) AiV(cj) 



where, for i = 1, 2 one has 
H,[t,N] 

Juj, L \ ""-^ ^ '-"^ ' 

duj. 

Hence, with the feedback laws 

u,{t,N):^-H,[t,N]yie{l,2}, (13) 

it follows that 
dC, 



dt 



[N{t)] = -ui{t, NY - U2it, NY ^ 0. (14) 



As in lBeauchard et ahl (|2011l ). we have the following re- 
sult. 

Proposition 2 For every initial condition Nq G 
H^((uJt:,UJ*), S^), the closed loop system (10), ( 13) has a 
unique solution N ([0, 00), ((w*, w 

*),R3)) such 

thatN{0) = iVo. 

2.5 Closed-loop simulations 

We assume here w* = 0, w* = 1 and we solve numerically 
the T-periodic system (1) with the feedback law {ui, U2) 
given by (2), (13). The closed-loop simulation is per- 
formed for t G [0, Tf], Tf = 20r and T = 2n/ {uj* - uj^). 
The w-profile 3 a; H- (x{t^Lo)^y{t,uj), z{t,uj)) is 

discretized {l,...,N + 1} 3 k ^ {xk{t),yk{t), Zk{t)) 
with a regular mesh of step ejv — " j^"* with N = 100. 
In other words, one has a set of discrete values {wi, i — 
1, . . . , -I- 1}, where w,; = + (i — l)ejv. 

We have checked that the closed-loop simulations are 
almost identical for N = 100 and N = 200. In the 
feedback law (16), the integral versus lj is computed as- 
suming that (x,y,z) and (x',y',z') are constant over 
](fc - \)eN,{k + \)eN[, their values being {xk,yk,Zk) 
and f '^+i-""-! VM^pj^ z,+,-z,^A ^j^g obtained 

differential system is of dimension 3(A^ +1)- It is in- 
tegrated via an explicit Euler scheme with a step size 
h = T/IOOO. We have tested that h = T/2000 yields al- 
most the same numerical solution a.t t — Tf — 20T. Af- 
ter each time-step the new values of (xk,yk, Zk) are nor- 
malized to remain in The initial w-profile Mq{uj) of 
(a;, y, z) G is given by xo = 0, yo = — \/ 1 — ' '^here 
zo = - cos(f ) + 0.05 (1 - cos(|) cos(cjf )). The desired 



final profile Mf{uj) is given hy Xf = 1 ~ Z/, Vf — 0, 
where Zf — — cos(Yg) -f- 0.1 (l — cos(^) sin(cj|-)) . 

The map R{uj) is constructed for the discrete set 
{wj,? = 1, . . . , A^ -I- 1}, in the following way. For 
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i — 1, one takes r^^uji) ~ Mo{uJi). Now choose a 
vector 9 among the vectors of the canonical basis 
in a way that {6, Mf{uji)) is the minimum value. 
Construct r2(wi) = jgj^^:^^ {6 A rsiuji)). Then 

one may take ri(aj2) = '''2{^2) A r^{uj2). Now, for 
i = 2, 3, . . . , iV + 1 one chooses ^3(0;^) — Mo{uji), 

e = -ri(w^_i), r2{uJr) = llflAraV.) ! ! ^ ^si^^i)) ^nd 

ri{uii) = r2{uJi) A r3(wi), and so on. The orthogonal 
matrix R{uj) formed by the column vectors ri, r2, r'a is 
then transposed to obtain R{uj). 

Figures 1 and 2 summarize the main convergence issues 
for these choices of initial profile Mq and of the desired 
final profile Mf. The convergence speed is rapid at the 
beginning and tends to decrease at the end. We start 
with £(0) w 0.1929. We get C{20T) « 0.0032. This nu- 
merically observed convergence is confirmed by Theo- 
rem 1 here below. 



Profile Stabilization 
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Fig. 1. Lyapunov function jC{t) defined by (12) and the 
closed-loop control (111,1*2) defined by (13) 

3 Main Result 

3. 1 Local stabilization 

Theorem 1 For every Mj G i/^((w*,a;*),§^) with 

(M/(w),e3) /0,Vwe (15) 

there existsSi > such that, for every Nq G if^((aj,, cj*), S^) 
with \\Nq -I- ealljyi < Si, the solution of the closed loop 
system (10), (13) with initial condition N {0 , uj) — No{uj) 
satisfies N(t) —63 weakly in H^{lo^,,lu*) when 

t -t-00. 

The above theorem has the following corollary. 

Corollary 1 For every Mf e 7Ji((a;*, w*), S^) with 
(15), there exists 82 > ^ such that, for every Mq S 




- M(T,) 

-M. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Fig. 2. Initial [t — 0) and final {t — Tf) u-profiles for x, y 
and z solutions of the closed-loop system 1 with the feedback 
(2), (13). 

i/i((w*,w*),§2) ^ith \\Mq - MfWni < S2, the so- 
lution of the system (1) with the initial condition 
M(0,w) = Moioj) and the feedback law given by (2), 
(13) satisfies M{{2kT)+) Mf weakly in when 
k — >■ -foo. In particular. 



hm ||M((2fcr) + , .) - M/|U..(„.,^.) = 0. 

/c— f-t-00 



The remaining part of this section is devoted to the proof 
of Theorem 1 



3.2 LaSalle invariant set 

The first step of our proof consists in checking that, lo- 
cally, the invariant set is reduced to {—63}. 

Propositions For every Mf G i7^((aj*, w*), 
with (15), there exists (5 > such that, for every 
No e iJi((w,,a;*),§2) mth \\No + e^W < S, the map 
t H- C[N(t)] is constant on [0, +00) if and only if 
No = -63. 

Proof: Let us assume that t ^ C[N(t)] is constant. 
Then mi = 7/2 = and N{t,u}) = No{u}) (see (14) and 
(10)). Thus, for every j € {1, 2} and t G [0, +00) 







/:.*[(^oH,(|£(i,^)e,)AiVo(^) 
+ (^63, (^F{t,Lu)ej'^ AiVoH 



(16) 



For t e [0,T], ait) = t so F{t,Lj) = Er=o '-^R^S" 
and ^{t,u) = E HrR'i'^)S' + E \^Riu:)SK 



k=0 



fe=l 
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Developing (16) in power series expansions of t and using 
(5), we obtain, for every j G {1, 2} and k ^ 1, 

+ (ea, (^i?(w)e,) A iVo(w)) = 0. 

By linearity, the following equality holds, for every Q G 
R[X] and j e {1,2} 



J (Kiu), {Q{u>)R'{w) + Q'{uj)R{co))e, A iVoH 



Q{uj)R{Lo)ej A No{u}) ) dco = 



+ (63, 

Thanks to the density of polynomial functions in 
w*), C), the previous equality holds for ev- 
ery Q G H^{{u}^,uj*),M.). Let us recall the relations 
{X,Y AZ) = {Y,Z AX) and {MX,Y) = {X,M^Y), 
VX, Y,Z and M e Al3(M), where M"^ denotes the 
transposed matrix of M. Then, the equality (17) may 
also we written 



/ {Q{ij)ej,R'{ioy[No{io) A iV^(a;)]) 

+ {Q'{uj)ej,R{ujy[No{uj)AN;,{uj)]) 
+ {Q{Lo)ej,R{ioy[Noioj) A es]) dio = 0. 



(18) 



for every Q G ff^ ((w* , w*), R) and j G {1, 2}. By linear- 
ity, we deduce that 

/ {Q{oj),R'{ujy[No{uj)ANl,{uj)]+R{ojy[No{oj)Aes]) 
+ (Q'iuj), R{uj)'^[No{uj) A iV^(a;)]) duj = 

(19) 

for every Q G i?^((w*, w*), V) where V := Span(ei,e2). 
Let P : ->■ V be the orthogonal projection on V. The 
previous equality is equivalent to 

j ¥R{ojy[No{uj) A es - Nq{oj) A N(^'{oj)] = in H'^ 
\ R{u}y[No{u}) A Ar^(w)] = at a; = and oj*. 

(20) 

Here, H ^ denotes the dual space of i?Q (w*, w*) for the 
L^-scalar product; the first equation has to be under- 
stood in the distribution sens. Thanks to ||A'o(w)|| = 1, 
we have ||i?(w)^[A^o(^^) A A^^(a;)]|| = l|A^o('^)ll- Thus, the 
second line of (20) is equivalent to Afg = at w* and 
CO*. Notice that Pi?(w)^|v is bijective on V for every 
u) G [a;*,a;*]. Indeed, thanks to (7) and (15), we have 



Range[Pi?(w) ■ 



: Ker[Pi?(a;)|v] 

: G V; i?(w) G Res} 

: VnMM/(w) = {0}. 



Moreover, (Pi?(w)"^|v)"^ G H^, thus (20) gives 

[ -A^^' A es + ^0 A es = (? in H~\{uj„uj*),Y), 



(21) 



[ A^Q A es = at lo»,uj*, 
where 

g{io) := -(Pii(a;)^|v)-'Pi?(a;)^K'(a;)A(Aro(a;) + es)]. 
Therefore, there exists Ci = Ci(w*, w*) > such that 

\\No{u)Aes\\Hi^C4g\\H-i. (22) 

Thanks to (8), there exists C2 = C2(w*,a;*, ||M/||jji) > 
such that 

\\g\\H-^^C2\\No + es\\j,^. 

When A^'o is close enough to — es in H^, then || A^o Aes||_f/i 
and ||A^o + es||//i are equivalent norms and then (22) 
gives 

||A^o(w) Aesllffi C3||A^o(t^) Aesll^i 

for some constant C3 = C3(w*,w*, ||M/||jji) > 0. This 
implies A^o A es = 0, i.e. A^o = ^63. □ 

Remark 1 ForMf = a, any constant function Nq with 
values in Span{e2,es) belongs to the invariant set (see 
(20)). Thus, an assumption of the type (15) is required 
for our strategy to work. 

3.3 Convergence proof 

For the proof of Theorem 1 , we need the following result . 

Proposition 4 Take Mf G i/i((w*,a;*),§2) and R G 
il^((a;*, w*), 50(3)) as in Proposition 1. Let {N'^)n£n a 
sequence of H^{{lo,,uj*),E^) and G iji((cj,, w*), S^) 
such that weakly in and N^^ Nj^ 

strongly in L^. Let a G [0, 2T] and (T„)„gK he a sequence 
o/[0, 2T) such thatTn — > ce. LetNn (resp. N^o) be the so- 
lutions of the closed loop system (10), (13) associated to 
the initial condition Nn{Tn) = (resp. Nao{a) = N^). 
Then, we have Nn{t) Noo(t) weakly in , Vt > a, 
and Uj[t, Nn{t)] -)■ Uj[t,N^{t)], > a,Vj G {1,2}. 

Proof: The sequence {N^)neN is bounded in and 
C[Nn{t)] sC C[N^], for every t G [r„,+oo) and n G N 
so there exists Mq > such that ||A^„(t)||^i ^ M.o, 
for every t G [t„,-|-oo) and n G N. The function t G 
M Fit,.) G i?i((u;,,a;*),A^s(K)) defined by (11) is 
continuous and 2T-periodic, thus, there exists A4i > 
such that ||F(t, .)||jji ^ Mi, for every t e R. Thanks 
to (13), we have \uj[t, Nn{t)]\ ^ 2MiMo, for every 
n G N and t G [t„,+oo). We deduce from (10) that 
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< iMiMo for every t G [r„,+cxD) and 

n G N . The end of the proof is as in IBeauchard et alJ 
(120 111) . □ 

Proo f of Theorem 1 : The proof is as in lBeauchard et all 
(j201l[ ). One may replace Barbalat's lemma by the 
Lebesgue reciprocal theorem, in the following way. 
Thanks to (14), t H> u^[t,N{t)] belongs to L'^{0,+oo), 
thus, for any diverging sequence of times /c„, the se- 
quence {t G (0, +oo) Ui[2k„T + i, N{2k„T + t)])nGN 
converges to zero in L^{0, +oo). Therefore, there exists 
a subset TV C (0, +00) with zero Lebesgue measure such 
that Uj[t, N{2knT + t) ^ for every t e (0, +00) - TV 
andjG{l,2}. □ 

4 Concluding remark 

Open-loop "impulse-train" control are combined with 
Lyapunov feedback to steer an initial profile 3 
u! Af (0, Lu) of the Bloch-sphere system (1) towards an 
arbitrary target profile [w:,,a;*] 3 uj t-^ Mf{uj). Conver- 
gence is proved to be local for any target profile belong- 
ing either to the south or to the north hemisphere. We 
guess that our convergence proof could be extended to 
the case where M f intersects transversely the equator 
and thus where Mf is not confined in only one hemi- 
sphere. 
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